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Abstract 

The pion decay constant and mass are computed at low temperature within Chiral 
Perturbation Theory to two loops. The effects of the breaking of Lorentz Symmetry 
by the thermal equilibrium state are discussed. The validity of the Gell-Mann Oakes 
Renner relation at finite temperature is examined. 
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1 Introduction 


The spontaneous breaking of chiral symmetry is believed to be a property of the 
strong interaction at zero temperature. If the temperature is finite different regimes 
appear. For sufficiently low temperatures chiral symmetry is still spontaneously 
broken, whereas for high temperatures it has to be restored according to asymptotic 
freedom. The way this restoration happens is not yet known and the transition 
temperature is estimated to be around 150 — 250 MeV (for a review of QCD at finite 
temperature see [jl|). 

The low temperature regime hadronic phase is dominated by the lightest particle 
occurring in the spectrum: the pion. The non-zero but small masses of the u and d 
quarks, which explicitly break chiral symmetry, make it a pseudo-Goldstone boson 
of the theory. Because of the lightness of these quarks, their masses can be treated 
as perturbations. 


When the system is heated, the hrst particles to be produced are the pions, 
whereas the states that remain massive in the chiral limit —>■ 0) are expo¬ 

nentially suppressed. The low-energy properties of the pions are essentially fixed by 
chiral symmetry. This leads to a wealth of low-energy theorems derived in current 
algebra at T = 0. For instance the Gell-Mann Oakes Renner (GOR) relation relates 
the pion decay constant and mass to the quark condensate and the quark mass (for 
simplicity all the quarks are put to the same mass m) |Q|: 


Ml 

m (0|gg|0) 


—1 -f 0{fh). 


( 1 . 1 ) 


A very efficient technique to analyze the corrections to these theorems is the use of 
effective Lagrangians. The low-energy effective theory of QGD is Ghiral Perturbation 
Theory (GhPT) [^ |^, p. At zero temperature and non-zero quark masses, the 
GOR relation ceases to hold Q. The goal of the present paper is to compute the 
temperature dependence of the quantities involved in the GOR relation and to see 
whether the latter still holds in the chiral limit. 


The thermodynamics of a hadron gas and the quark condensate below the chiral 
phase transition have already been studied up to three loops in GhPT 0. The 
question of the pion propagation at hnite temperature has been addressed: the self¬ 
energy is known to two loops in GhPT [0, but the pion decay constant to one loop 
only i. 

The thermal equilibrium state is not Lorentz invariant. In GhPT this first ex¬ 
plicitly shows up at the two-loop order. It is then important to compute up to 
this order to see what happens in the non-zero temperature case. For instance the 
appearence of two distinct pion decay constants — one for time and one for space — 
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begins at this order. This was already noted in ||^ and is a known featnre of some 
non-relativistic system snch as the antiferromagnet [pil|] . 

A short presentation of ChPT at hnite temperature in the Real Time Formalism 
is the subject of the next section.The two-point axial Green’s function is computed 
in Section 3, the pion decay constants and mass in Section 4. The GOR relation 
at hnite temperature is the subject of Section 5. It is derived from two Ward 
Identities relating the axial and pseudoscalar two-point Green’s functions. Some of 
the consequences of these identities are also examined. The regularization and the 
chiral limit of T-dependent functions appearing in the corpus are appended. 


2 Cool Chiral Perturbation Theory 


A hnite temperature ehective theory rests on the zero temperature one. In the path 
integral formalism going from one to another is essentially a change of manifold over 
which one integrates: a torus replaces a plane |^. The Real Time Formalism is 
used because we are investigating time-dependent Green’s functions. 


To make the presentation simpler, we will only give those results we consider 
necessary here. The reader is invited to consult ^ for a more detailed discussion 
on GhPT at T = 0. We will restrict ourself to the case where all quarks have the 
same mass m. 


GhPT is an ehective theory describing QGD at low energies. The iV-havour 
massless-quark QGD Lagrangian is symmetric under SU{N)ji x SU{N)l, the chiral 
group. It is assumed that a spontaneous chiral symmetry breakdown occurs, 

SU{N)r X SU{N)l ^ SU{N)v, 

whose Goldstone bosons are identihed as the pions. 

The QGD Lagrangian can be approximated at a given order in the momentum 
using an ehective Lagrangian expressed in terms of a held U G SU{N) which trans¬ 
forms linearly under SU{N)ji x SU{N)l, 

U ^ gnUgt, 

and contains the helds of the pseudoscalar Goldstone bosons, 

U = (2.1) 

where r“ are the generators of SU{N) and F is the pion decay constant in the chiral 
limit: = F{1 -|- 0{m)). 
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Coupling U with external fields and expanding the effective Lagrangian in powers 
of the external momenta and of quark masses, gives 

/leff = + • • • (2.2) 

The promotion of the global chiral symmetry to a local one requires the introduc¬ 
tion of a derivative which is covariant with respect to the external axial and 
pseudoscalar gauge fields. 

To get the desired 0(p®) accuracy, the tree, one- and two-loop diagrams of 
the tree and one-loop graphs of and the trees of are needed. 

We will restrict ourselves to the two flavour case {N = 2) and use the tradional 
notation of [^]: y = 2Brh I -|- 2iBp, = d^U — f{a^, U} and = Fd^ay =F 
d'^a^ — i[a^,a^], where B is proportional to the quark condensate in the massless 
quark limit. It is found that the parts of the Lagrangian that contribute to the axial 
and pseudoscalar two-point Green’s functions to 0(p®) are: 

1) Lowest order: 

+ yf/t + ytf/) (2.3) 

2) Second order: 

+L,{xU^ + ytl7)(V^17tV^t/) + L,{xU^ + ytf/)^ (2.4) 

+L^{XU^XU^ + X^Ux^U) + Lio(f/^Frf/T’LM.) 

+ H,{xh) 


3) Third order: the singular part of the two-point function we are interested in 
only receive T-independent corrections from the trees of £efr (see Section 3 for an 
illustration of this property at one loop). Hence won’t be explicitly needed here. 

When the temperature j3~^ is non-zero, the fields U{x) map the /3-dependent 
torus into SU{N). The generating functional of the connected Green’s functions is: 


exp 


{^iZ[a^,p ]}=J |dC/| exp |i y d^x drCeff 


(2.5) 


The integration extends over IR^ and a contour C. In the Real Time Formalism 
of quantum field theory at finite temperature, which has been thoroughly studied 
in 1^, one can choose different integration paths in the complex f-plane. We take 
the so-called Keldysh path shown in Figure 1. The functional integral extends 
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over field configurations with the boundary condition U{—oo,x) = U{—oo — il3,x), 
which is the familiar periodicity condition. A generalization of time-ordering has to 
be introduced: T is replaced by Tc, an operator which orders operators according 
to the occurence of their time-arguments on the contour. A Heaviside, 6c, and a 


Dirac, 6c, distributions are dehned on the contour |11 


The thermal propagator is a Green’s function on the contour: 

(Dc + M^) D/ 3 {x- y) = 6c{x-y). (2.6) 

The unique solution satisfying the Kubo-Martin-Schwinger boundary condition is 


Dy{T - T']UJk) = 


2uJk 






Oc{r' -t)Y 


(2.7) 


+ 


iWk{T-T') _|_ ^-j3uJk-iu}k{T-T') 


where nsix) = — 1) is the Bose distribution and ujk = ^k‘^ + M^. 


When ( ^.5| ) is used to compute Green’s functions with arguments on the real 
f-axis, only a part of Z[a^,p] is relevant: the Ci and C 2 contour segments are the 
important ones. 


The introduction of the usual 2x2 formalism allows us to rewrite the generating 
functional as (j = a^,p): 

Z[ji,j2\ = l[dUi][dU2] exp|* j ^UmD^^^^Un-V[U,] + V[U2]+jnUnY (2.8) 
where m,n = 1, 2, and ji{x) = j{x), j 2 {x) are independent sources. 


The real-time Green’s functions are generated by differentiating (|2.8|) with re¬ 
spect to ji and setting both ji and j 2 to zero. For instance the real time two-point 
axial Green’s function is given by: 


i{TAl{x)Al{y))T = 


Tre-^^ 


Tr [e-P^TAl{x)A 


5^ 


6a1^{x) 6a\^{y) 


Z[ai^, a2fj,,Pi,P2] 


. (2.9) 


ai(j=a2/i=Pl=P2=0 


In momentum space the propagator D™" reads: 


iDf{k) 

tDfik) 


tDfik) 

iDf{k) e- 


= A{k) + 2 ti 6{k‘^ — M^) n^dfeol) 

6 = 271 6{k‘^ — M^) nB{\ko\) 


Af3{k) 


( 2 . 10 ) 
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where A(fc) = ij + ie) is the usual Feynman propagator. 

The topological and combinatorical structures of the Feynman diagrams do not 
change, but the propagator has now acquired a matrix structure and there are two 
kinds of fields. In a real-time Green’s function only the type 1 field can appear on 
an external leg and the type 2 field plays the role of a ghost field. There are also two 
kinds of vertices, which are equivalent up to a sign. The two different fields interact 
through the off-diagonal terms of the propagator. 


3 The axial two-point Green’s function 


The hnite temperature two-point Green’s function of the axial current 

= 9(2:) y7(.T5?(i) 


(3.1) 


describes the dynamics of a pion in a strongly interacting gas (r“ are the SU{2) 
generators). Let — y,T) S°‘^ := i{T A“(|/)yl[)(|/))T. In general the two-point 

axial Green’s function at finite temperature can be written as: 




T) + q^,(3^{q, T) + quP^iq, T) + T) 

ql-^{q,T) 


^p%{q^T) (3.2) 


It contains a singular piece and a finite part, p^^{q, T). To get this Green’s function 
to order in S'f/(2)-GhPT (two flavours), one needs to compute the Feynman 
graphs shown in Figure 2 with the effective Lagrangian ( p.2|) . 


The singular parts of the tree diagrams are temperature independent. They just 
contain the zero temperature Feynman propagator A(g). Most of them have already 
been computed in Q. The other components of the matrix propagator generate 
imaginary parts of the finite term of the two-point Green’s function. Thus they 
are not interesting in our context. They will only be mentioned in the one-loop 
computation given below and ignored elsewhere. 


The T-dependence of the one-loop graphs is rather simple: only properties of 
the thermal propagator at the origin appear in that case Q. 

At lowest order the pion mass is 


:= 2rhB 

(3.3) 

and the various O(p^) graphs are: 


+ 16M% - 2A^(0)) 

(3.4) 
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(3.5) 


G2;«)(9,r) 


= I 




+ ['?^ Qfj.i' ~ QfiQv] — 2Lio) 

/3(0) + 3 

4g2 - m2 


QiiQu ( ^A^(O) + 32 L 4 M 2 ) Ap{q) 
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A/3(0) 


(3.6) 


+M2(-16L4g2 + 32LeM^ + IGLgM^) (A^(g))2 - {tDf{q)) 


The relation (□ + M‘^)Af^{x) = i5^^\x) has been used to simplify the expressions. 


At hrst sight, the term in square brackets in (|3.6|) seems to be ill dehned: it 
contains products of Dirac distributions with the same argument. But the equation 
of motion implies that 


dM^Afsix) = i^Dfix) 


dM‘^ ^ 

i j d’^y [Df{x - y) Df{y) - Df{x - y) Dj 


(3.7) 


Hence the products of thermal statistical weights at intermediate stages of the com¬ 
putation disappear when the different types of vertices are combined. This is a 


generic feature of the Real Time Formalism [11 


Moreover the terms coming from the thermal parts of Ay{q) and the alike only 
contribute to p^^{q,T), the hnite term of the two-point function. To get the pole 
and residue of the latter, it suffice to replace A^(g) by the corresponding Feynman 
propagator everywhere. This is the case to all orders in the perturbation. 


The functions appearing in (|3.2|) at first non-leading order are already known 


f m2 

7") = q"^ + m 2 j 1 + [—I 6 L 4 -|- 32L6 -|- IGLg] 

m2 

+ ^MO) + O{p^), 


(3.8) 


a^{q,T) = 
= 

ptM^T) = 


2F2 

M2 


1 + ^ [I 6 L 4 : 


-2A^(0) + O(p®), 

0(p®), 7A = 0(p®), 

(f2 + 16M2L4 - 2A^(0)) 

+ [q^ 9tiu - Qi^Qu] (4Mi - 2Lio) 0{p^). 


(3.9) 

(3.10) 


Only p^^{q,T) really depends on the momentum. As already said, the Lorentz 
symmetry is not explicitly broken at this order. 
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Because A; 3 ( 0 ) contains the Feynman propagator at the origin, it diverges. As 
usual in ChPT, because it preserves the symmetry of the theory, dimensional regu¬ 
larization is used [§, |[: 

A^O) = / 0jsMk) = 2MU + Um(T), (3.11) 

where 

^ ^ - illn to + r'(l) + 1] + ln(") + 0(d - 4)) . (3.12) 

and 

r 

Mm{T) = J - m 2 ) hbM. ( 3 . 13 ) 


The renormalization of the theory at hnite temperature has to be the same as 
at T = 0. This is obviously here the case: no divergence occur in the temperature 
dependent part of the two-point function. The scale independent parameters used 
to renormalize the theory are: 

Li = Li-^iX, Hi = Hi-6iX, (3.14) 


where 71 = 1/12, 72 = 1/6, 74 = 1/4, 70 = 3/32, 73 = 0, 710 = 2/3, = 1/3 and 

^2 = 0. In ChPT at T = 0, the expansion of the pion mass and decay constant at 
O(p^) can be expressed in terms of F, M and Lf. 


Ml 

= M^ 

F2 

TT 

= F2 


M2 r 

■ JF _ 
M 2 - 
Y2 L 


+ 0 (/), 


16L4 




(3.15) 

(3.16) 


The 0{p^) graphs are a bit more complicated than those of 0(p‘^). However, 
only the T-dependent terms are of interest here. Thus the diagrams (6.23 — 28) will 
not be explicitly given 


see 


12, 13, 14] for an analysis of zero temperature ChPT to 


0 (/)). 


Some of the two-loop graphs, (6.1 — 3, 6.13 —16), are just products of lower order 
ones. 


The diagrams (6.4 — 6 , 6.17 — 19) involve essentially the same elements as those 
occuring at one loop. Only their vertices are different: they contain some low-energy 
coupling constants and derivatives. For instance the graphs (6.5) and (6.18) give: 

8 ? 

(3-17) 
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~ 1 (S? ~ r rl^lc ~ r 

j ^/ 3 (^) {( 2^1 + 4 ^ 2 ) kq {k^q^ + k^q^) ( 3 . 18 ) 

-\-[Qk^Li + 2k^L 2 — —M'^L^) q^iqu^- 


The graphs (6.7 — 9, 6.20 — 22) are special ones for the Real Time Formalism. 
They are the only diagrams in the whole set which contribute to the pole and residue 
of the two-point function and contain type 2 helds. Their role is very important for 
the consistency of the theory [^. The graphs (6.8) and (6.21) are taken as examples: 


(^( 6 . 8 ) 
^ flL/ 


4 


q^q^ A(g) A^(0) 



{iD'phif - {iDf(k))‘ 


(3.19) 


where the fact that the type-2 tadpole is the same as the type-1 tadpole has been 
used, and 


(^( 6 . 21 ) 


128M2 ~ ^ ^ 

3^2 


r d‘^k 

J {27iY 


{k‘^U - 2 M^Lg - M^Ls) 

{iDf{k)y - (iDf{k)y 


(3.20) 


As a consequence of ( p.7D , the integrals in ( p.l9| , p.20D are essentially 9 m2 A^(0). 
In d = 4, this quantity contains a singular piece: 


^M2A,g(0) — 2A -I- dM^AfiMiT). 


(3.21) 


Finally, we turn to the genuine two-loop graphs. These are the closed-eye (6.10), 
the guimbard (6.11) and the sunset (6.12), given by: 


{(3^1^ + 3/ c 2 /i ~ 2 g^) {3kii, -|- 3k2u — 2qY} , 

^ q^,H(l) [ ^ ^ A^(fci) A^(fc 2 ) A^(g - h - k2 


(3.22) 


g^2 J ('27r)2'i 

i2qy — 3/cii/ — 3k2i/) -|- + 4:kik2 -|- 2 (fci -|- /C 2 )? — 2 q'^| 

/i <—^ u ), (3.23) 




?/2gi^A2(g) J ^ Afs{ki) Afs{k 2 ) Afs{q - ki - /C2) 

3M^ + (fc? + kl + Akik2 + M^ + 2{ki + k2)q - ^q^^ . (3.24) 

























One readily sees that the polynomials appearing in the integrands are essentially 
four-point functions at tree level. For instance the square of the tree-level isospin 
averaged vr-vr scattering amplitude appears in the sunset. These integrals may hrst 
appear complicated, but using the symmetry properties of the integrands they can 
be expressed in terms of two independent integrals (cf Appendix A for more details): 

G!.y“'(9.r) = ~ I(q,T) + T V(«,T), (3.25) 


(3.26) 

-Aq^qPI^piq, T) + Aq^q^A^ 


= -M.^{(5AC-8«‘)/(«,r) 


18F2 


(3.27) 


+72gW(g,T) - 12(M2 - 3g2)Aj(0)} . 


The functions I{q,T) and Ip_^{q,T) are dehned as: 

Hq.T) = »/ A^ik) A„(fc;) A^(q -k,- h), (3.28) 

^ ^p{q -ki- k2) kip ki^. (3.29) 

These integrals have to be regularized because they are divergent in d = 4. Their 
hnite parts are determined by the four functions A/m(T), Np^{M,T), I{q,T) and 
Ip^{q,T), given in Appendix A. 


The end result for the various terms appearing in the representation (|3.2|) reads: 


A/f 2 

■(g,T) = q^ + 


2 Ff 


1 


Mi 


ATmAT) + -fATMAT) Om^MmAT) 

11 - 

— ^UAAT) + /(g,T) - gV«..(g,T) 


(3.30) 

+ 0 (p«), 


a 


'(q,T) = F;-2Mm,(.T) 


(3.31) 


F 2 t 


M^ UmAT) - MIMmAT) OmiMmAT) + 2Ai2 (T) + 0(/), 
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/5^(g,T) 


(3.32) 


= i + 0{p^), 

■^TT 

7i = 0{p^% (3.33) 


P^,M.T) = R^^{q)+g^^a {q,T) 
1 


F 2 


K 




(9. T) + ( 


qi,q,y - 9t,^ q 



+ 0(p«). 


(3.34) 


The function 

R^,y{q, T) := L T) + 44^(g, T), (3.35) 

and different combinations of the renormalized low-energy coupling constants were 
introduced to lighten the expressions: 


( L — 32Li^-|- 641/2 — Yg^ 

< Is = -48Li - I 6 Z 2 + 48 Z 4 - sole - iOLg + ^ (3.36) 

[ La = 48Li -|- I6L2 — 24L4 — . 

Expressed in terms of the coupling constants dehned in [Q, these are given by 
ri = (ri + 4r2-^)/i27r2 

Lh = (-24/1 - 16/2 +15/3 + 12/4 + f)/1927r2 (3.37) 

[ La = (6l + 4/1 - 9U - |)/487r2. 

Note that a^{q,T) does not depend on the momentum at this order. 


An integral representation can be given for the different functions of temperature 
and momentum appearing in the previous expressions. AV(7") has already been 
introduced in (|3.13| ). The other function related to the properties of the propagator 
at the origin is T). It is given in Appendix A together with the regular parts 

of the momentum dependent functions I{q,T) and I^^{q,T) dehned in (p.28|,|3.29D . 
All these functions depend on the ratio M/T in a non-trivial way. 


The expressions ( |3.3CI| - |375^ ) contain all the contributions to the hnite tempera¬ 
ture axial two-point Green’s function to 0{p^). Both M and T count as quantities 
of 0{p). The dependence of the functions H(g,T), a^{q,T) and P^{q,T) on the 
external momentum begins at this order. There is an important qualitative change 
between the one-loop and the two-loop results: the way the functions involved de¬ 
pend on go and on q are now different, rehecting the breaking of Lorentz symmetry 
by the heat bath. As it must be, the T 7 ^ 0 renormalization is the same as the one 
at zero temperature |^. 
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4 The pion decay constants and mass at finite T 


From the two-point axial Green’s function some interesting quantities can be derived: 
the pion decay constants and mass. In the T = 0 case the self-energy is dehned as 
the pole of G^u{(l,T = 0) in the go complex plane and the pion decay constant as 
its residue at the pole position. 

The dispersion curve determines the position of the pole in the go"Plane: 


go = G(g,T). 


(4.1) 


In our case a non-trivial momentum dependence occurs at 0(p®). Therefore the pole 
position at the same order can be obtained by replacing go by Uq = that 

is: 


n\q,T) 


{qo,q,T) _ +0(/). 


qo=ujq 


(4.2) 


A possible dehnition of the mass is := Re fl{q,T) ^ i.e. the real part 

of the pole. With the Ls dehned in ( p.36|) , the ChPT result reads: 


M 2 (T) = i + 


1 


Ml [UmAT) OmiMmAT) 


Ml 

6 


11 


Re[I{q,T)]--Afi,^{T) 


(4.3) 


-LiVoo(M^,T)-4Re[Joo(g,T)] 


+ 0(/). 


qo=M,q=0 


This result agrees with 0, where a somewhat different representation is used. 
The imaginary part of the pole, which determines the damping rate of the pions, 
has been thoroughly studied in the same context: its mean approximately behaves 
like ryPl above 100 MeV §. 


To extract the residue at 0(p®), S(g,T) has to be expanded around the pole 
position: 

.,dE{q,T) 


^(go, pT) = n (g,T) + (go - ujg) 


dqo 


qo=uJq 


+ ... 


(4.4) 


The thermal equilibrium state is invarian under spatial rotations, it implies that 


Koi(g,T) = qi K{q,T), 


(4.5) 


11 















(4.6) 


thus P^{q,T) can be rewritten as: 


(3^{q,T) = go Pf{q,T) and (3f{q,T) = q, Pf{q,T). 
The axial two-point function then takes the form 


where <hj^j,(g,T) is a hnite term and 


qi-n^{q,T) ’ 


/o(g,T) = qoF,{q,T) 

Mq,T) = q.FMT). 


(4.7) 


(4.8) 


The ’’temporal” residue at the pole position is given by 

1 dE{q,T) 


F,\q, T) = (a^(g, T) + 2pf{q, T)) (l + 


2ug dqo 


qo=uJq 


+ 0{p^ 


(4.9) 


The difference between Pt{q,T) and Ps{q,T) makes the difference in Ft{q,T) and 
Fs{q,T). 


The pion decay constants can naturally be dehned as F^’^(T) := Fs^t{q,T) 
The ChPT computation gives 


q=0 




Fi 


Af- {MmJt) s„w„,(r)) 


/Vf3 d _ 


+ 4 / 00 ( 5 ', F) — 2Mt,— —/oo(g, T) 
oqo 


+ 0{p^ 


qo=M,q=0 ^ 


and the difference between the two pion decay constants is 


F^(T)-F;(T) 


1 


i L (4Af„„(Af„T) - APMmST)) + \K,m (4.11) 


4 - 16 - 

--M2j(g, T) + yJoo(g, T) - M.k(g, T) 


+ 0(p®). 


qo=M,q=0 


The various functions of temperature involved in the expressions above are given in 
Appendix A. 
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The T-dependence and the evolution along the perturbation expansion of the 
mass and the ’’temporal” pion decay constant are displayed in Fig.3-4 in the case 
of the physical pion mass and pion decay constant: ~ 140 MeV and ~ 93 

MeV. For both and F^{T) the third order corrections have the opposite sign 

as those of the second order. At one loop, the mass is enhanced by the effects 
of the temperature, but the two-loop corrections bring it down (Fig.3). This may 
be a reflection of the fact that at T = 0, the hrst correction is negative for M 
(|3.16|) . Exactly the opposite happens to the pion decay constant (Fig.4 and (|3.16|) ). 
The difference between the ’’temporal” and ’’spatial” pion decay constants is always 
positive (cf Fig. 7 in the next Section). 


5 Massless quarks 


In our problem, the small and hxed T case is equivalent to the limit T 3> Mt^\ 
only the ratio M^/T is relevant. In the chiral limit, that is when the quark masses 
tends to zero, tends to zero. The expressions of the pion decay constants and 
mass are much more readable (F ~ 88 MeV is the pion decay constant in the chiral 
limit 0 ): 


Ml{T) 


Ml 


1 + + 


m=0 


24F2 


36F4 

1927r2 


19 , T 

480 /i 

{L\{p) + 2L'M) 


(5.1) 


+ 0 (T“), 


Re(F^(T))' 


FI 


= 1 


2^4 

+ 


m=0 


6F2 36F4 

1927r2 


7 T 

-F - In - 

60 p 


(5.2) 


+ 


(Ll(/i) + 2L”(p)) 


+ 0(T®), 


{Fi(T) - f;(t)) 


2 

F^ 

27F4 

m=0 

15 


1927r" 


(LK/i) + 2L”(/i)) 


(5.3) 


+ 0 (T®). 
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Where the number 


1 45 

K = \n2 + -r'(l) + ^C'(4) - 1.05 ~ -0.68 (5.4) 

2 TT^ 

contains the Euler gamma and Riemann zeta functions and a contribution from the 
integrals that had to be numerically evaluated, /i is the regularization scale used 
in (|3.12|) and E[(/r) = Lj + 7 * ln^/( 47 r)^ are the scale-dependent renormalized 
effective coupling constants. The scale dependence and the logarithmic divergences 
of the individual terms appearing in (b.l[-b.3|) cancel, as they have to. 


Our expressions for the pion mass and decay constants can always be written in 
the form 


M 2 (T) 


M 2 


= 1 + 


m=0 

2 


_ 

24 F2 T + ^ 


F3 


= 1 


rp2 rp4 A 

-s +-J In — + on*), 

6F2 3eF‘ T ’ 


(5.5) 


(5.6) 


m=0 


Re{Fl{T) - F^{T)) 


F^ 


^^ln^ + 0(T®). 

27 pi T ' 


(5.7) 


m=0 


Where AM,t,A are various scales which sizes are determined by the numbers and the 
values of the coupling constants appearing in ( |5.1| - |^ ). We take the recent two-loop 
evaluation as reference: h = —1.7 and I 2 = 5.4, i.e. 327r^(Li-|- 2 L 2 ) — 1.66. The 
scale we hnd are rather big compared to the ones usually involved in ChPT 


Am 

~ 1.9 GeV 

(5.8) 

At 

~ 2.3 GeV 

(5.9) 

Aa 

~ 1.8 GeV. 

(5.10) 


In the range allowed for the temperature. In is positive. As already seen 

in the physical pion mass case, the third order corrections have the opposite signs 
to the second ones in M^(T) and F^{T), whereas the latter is bigger than F^{T). 
This is shown in the Fig.5-7 both in the chiral limit and in the physical case. 


To see what becomes the GOR relation at order T^, we need the quark condensate 
to 0{p^). It has been computed to 0(p®) in |^. In the general case it reads: 


m{qq)T = m{0\qq\0) +-M^AfM^{T) 


Ml 

F 2 


- 24 T 4 -|- 48Tg -|- 48T8 -|- 


647r2 


MiNu.m 

a/Lf2 Q 


(5.11) 


+ 0 (/), 
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thus in the chiral limit 




m=0 


= (ote|o>{i-^-^) + o(r«). 


(5.12) 


Which means that a modihed GOR relation can be written down for massless quarks 
at hnite temperature: 


lim ^h(I)MthI))!l = _i + 0(T«) 

m^o m{qq)T 


(6.13) 


This was already noticed in 0 and and is a consequence of the Goldstone 
theorem at hnite temperature as will be seen in the next Section. 

As expected, the essential characteristics of the low temperature behaviour of 
the three quantities examined are already present when the quark masses are sent 
to zero. Moreover, in this approximation a nice interpretation of the difference 
between the two pion decay constants and the imaginary parts of the residues of the 
two-point functions can be given. The next Section will come to that point. 

Now, the symmetry groups occuring in the 0(4) Linear Sigma Model (LSM) 
are the same as those of QGD with two massless havours. Hence the effective held 
theories of these two systems are identical, only their ehective coupling constants 
dihers. This implies that the results obtained above are valid as they stand also for 
the 0(4) model. In particular, the temperature expression of the pion mass contains 
a specihc logarithmic contribution at order T^. As it is absent in the LSM calculation 
described in |^, their result is not complete at order T^. The coupling constants 
LKfi), L 2 {fi) can be evaluated for the LSM with [H. The logarithmic contributions 
may be viewed as arising from a temperature dependent ehective coupling constant 
of the LSM. Indeed, using the formula (34) in P 


M^AT) 

Ml 


= 1 + 


y 2 




6 F2 15 F2 


m 


2 ’ 


(5.14) 


and replacing rrio- with ma(T) dehned by 


m2(T) 




247r2 F2 


In 


T 


a rrin 


(5.15) 


one recovers our result, provided that a ~ 0 . 68 . 

The same reasoning may be applied to the expressions found in P for both F^’*(T). 
Again the logarithms occuring at order are missing. These expressions are com¬ 
patible with GhPT if is replaced by m„{T) like in (^.14 , ^T5 ) but with a ~ 0.62. 
In the notation used above the diherence between the two logarithmic scales arises 
from Am 7 ^ A^. This diherence also manifests itself in the temperature dependence 
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of the quark condensate at order T^. If the formulae in were correct, the tem¬ 
perature expression of the quark condensate would not contain a term of order T*^, 
in contradiction with the old calculation described in |p, where the result is given 
up to and including contributions of order T® In T. 


6 Gell-Mann Oakes Renner relation at finite T 


To see why the Gell-Mann Oakes Renner relation has to take the form ( b.l3| ) when 
T 7 ^ 0 and to compute the corrections in the quark mass is the main purpose of this 
Section. A possible way to do this is to go back where it originates at T = 0. It 
may be derived from two Ward Identities involving the quark condensate and both 
the axial and pseudoscalar two-point Green’s functions. As already mentioned the 
difference between the T = 0 and the T ^ 0 cases is a change of manifold in the 
path integral formalism. Hence the derivation of Ward Identities from the generating 
functional of QGD at hnite temperature go through the same steps as at T = 0. 
In the hnal equalities the vacuum expectation values of the involved operators are 
just replaced by their thermal average. These Ward Identities will be used both 
as a consistency check of the whole computation and as a source of the generalized 
relation. They will have byproducts which are going to clarify the meaning of our 
results. 


To construct the hnite temperature sisters of the Ward Identities which lead to 


the GOR relation at T = 0, one has to use the axial-vector current ( p.l|) and the 
pseudoscalar density 

P“(a;) = g(x)yZ75g(x). (6.1) 

With the QGD Lagrangian, the following Ward Identities are obtained: 

= 2mG'^(g,T), 


qf^G^{q,T) = 2mG{q,T) + - {qq)T, 


( 6 . 2 ) 

(6.3) 


where G^y{x,T) is the axial two-point Green’s function dehned in Section 3, 
G^{x - y,T) 5“'’ := i{T ARa:)P^(|/))r and G{x - y,T) 5“'’ := i{T P^{x)P\y))T. 


Note that at T = 0, both identities are fulhlled by GhPT. The hrst one implies 


that 


27hG, = M‘F,, 


(6.4) 


and together with (|6.3|) in the chiral limit, one hnds the GOR relation at zero 
temperature. 
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The new two-point functions are represented in the same way as the axial one 

(E3): 


G,(g,T) = d>;,^(g,T) 


G(g,T) = $^(g,T)- 


f,{q,T)g{q,T) 
qi-n^{q,T) ’ 

9^q,T) 

qi-n^{q,Ty 


(6.5) 


( 6 . 6 ) 


r2(g, T) is dehned as in ( |4.2|) and the residue g(q, T) similarly as the one of the axial 
two-point function ( [4.9| ). Again the strength of the coupling of the pseudoscalar 
density to the pion is dehned to be Gt,{T) := g{q,T) ^ like in the axial case. 


The hrst Ward Identity implies that 

n\q,T) fyq,T) - f f^T) = 2m g{q,T). 
At g = 0, it generates a relation very similar to (|6.4|): 


(6.7) 


2rhGAT) = ff{q,T) 


( 6 , 8 ) 


Then the second Ward Identity together with ( |6.7|) implies that 

- F,(g,T) g(g,T) (g,T) = 2m$^(g,T) + (6.9) 

The quark condensate is independent of g. the previous equality can thus be evalu¬ 
ated in the chiral limit and at g = 0 (or conversely). Together with (|6.7|) it gives a 
generalization of the GOR relation at hnite temperature in a form very close to the 
one at T = 0 (and to all orders in the perturbation theory): 


niq,T) _ {Fmr 

lim --= -l. 

rh{qq)T 


( 6 . 10 ) 


Because the quark condensate is a real quantity, the following relation must hold: 

\2 ■ 

- = 0 . ( 6 . 11 ) 


Im 


lim 

m—5-0 




m 


Our result explicitly verihes this property at 0(p®). 


In order to compute the hrst corrections in the quark mass to the GOR relation at 
hnite temperature, it is enough to compute the pseudoscalar two-point function. It 
requires the same building blocks as those of Section 3. Therefore only the end result 
will be given here. The graphs involved are the same as those appearing in Figure 2. 
The vertices are in general diherent and the diagrams (2.1,4.1,6.4,6.7,6.20) are 
zero, because is only linear in the pseudoscalar external held. 
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To present the result a representation similar as the one we used for the axial 
case in Section 3 is used: 


G{q,T) = 


a^{q,T) 


+ P^(g,T) 


%^-S(g,T) 

The various functions appearing in the previous expression are found to be: 


( 6 . 12 ) 


a 




1 

'f2 


Ml Ip Mu.(T) - -MINm,(T) 




+ Oi/), 


(6.13) 


p’’{q,T) = R‘’(q) + -^ 


Lp UmJT) --UqX) 


+ 0(p«), 


where two new combinations of coupling constants have been used 


Lp = -48Li - I 6 L 2 + 104L4 - 224L6 - 112L8 + 

Lp = 32L6 — SLg — 3 ^. 
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(6.14) 


(6.15) 


The expression for S(g, T) is the same as ( p.30| ). This was a hrst test for the whole 
computation: the poles of the two Green’s functions under consideration must be 
identical. 


The corrections to the massless quark world can be computed. The hrst one that 
appears is 


G(g,T) (F^(T)) ^ ^ _ 

--= - 1 + 2 ^ ^^oo(g,T) 

m{qq)T dqo 


+0{M^ln^,p^), (6.16) 


qo=M,q=0 

the expansion of the last integral in terms of the pion mass is given in the Ap¬ 
pendix B. We get: 

\2 ■ 


Re 


Q{q,T) {F^{T)y 


m{qq)r 


(6.17) 


A small deviation linear in the pion mass is the hrst that appears. It is of the 
order of 6% at T = 100 MeV and around 20% at 150 MeV. This is in contradiction 
to the QCD sum rule result obtained in [0, where the hrst correction is quadratic 
in ikL. 


The Ward Identities (|6.2| , |6.3|) can be used as a consistency check of the whole 
calcualtion. The expressions for the quark condensate (|5.11|) , the axial (|3.2|) and 
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the pseudoscalar ( |6.12D two-point Green’s functions must fulfill the identities. It is 
here the case, irrespective of the specihc forms of the different independent functions 
involved in them. This is also the case if one includes the imaginary parts of the 
hnite terms that were ignored (see Section 3). 


The relation (|6.7|) is very instructive in the chiral limit. It implies that 


n{q,T) 


e^(g,T) 


(6.18) 


Thus the speed of the pions in the chiral limit is 

+ 0(/). (6.19) 

7h=0 

It has to be smaller than the speed of light. Hence the difference between the real 
part of the two pion decay constants has to be positive (this remark was already 
made in 0). Our result verihes this property. But looking at Fig. 7 or at (^.7|) , 
one sees that the square of the speed tends to zero at T ~ 160 MeV and even turns 
negative beyond that point. This is of course not allowed and it dehnes a natural 
limitation of the 0{p^) ChPT expansion. 


m=0 


= 1 - 




F^ 


Finally, because all the quantities involved in ( 6.18 ) are complex, this identity is 
in fact a system of two equalities containing six unknowns: the real and imaginary 
parts of the pole and of the two residues. Hence in the chiral limit, lmF^’^{T) are 
completely determined by the other four quantities. Three of them are given in this 
article, whereas the imaginary part of the pole has been deeply studied in [0. As 
a consequence, the physical content of the imaginary part of the residues encodes 
that of the other quantities involved, which physical interpretation is well known. 


7 Summary and Conclusion 


The dynamic of a pion travelling through a gas of pions at low temperature (up 
to circa 150 MeV) can be computed with the help of the ChPT Lagrangian. The 
knowledge of the effective mass is important to understand how this happens. The 
two different effective pion decay constants that appear because of the breaking of 
Lorentz symmetry by the equilibrium state are also quite meaningfull. They were 
derived here to a accuracy performing a two-loop calculation in ChPT. The 
temperature dependence of these interesting quantities is small. This is due to the 
size of the pion decay constant at T = 0 which governs the ChPT expansion. 

When the quark masses are sent to zero, contrary to the result obtained within 
the Linear Sigma Model in ^j, we hnd a logarithmic dependence in the temperature 
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of the three mentionned observables at the two-loop order. The difference between 
the ’’temporal” and ’’spatial” pion decay constants is positive. This is related to 
the fact that the velocity of the pions in the gas is smaller than the speed of light. 
The Gell-Mann Oakes Renner relation is still satished, which is a reflection of the 
Goldstone theorem at hnite T. 


For realistic qnark masses, the temperature dependence is quite obscure: it is 
contained in complicated integrals and some pictures are needed to see what hap¬ 
pens. The behaviour of the pion decay constants and mass for the physical pion 
mass is not very different from the massless case. 


The hrst corrections to the GOR relation are linear in the pion mass, they have 
a small magnitude. This does not agree with the result given in [T^ using QGD 
sum rules. 


Finally the presence of massive particles in the thermal equilibrium state deserves 
a comment. Their effects have been carefully analysed in P]. The lightest ones that 
appear in our case are iF(500) and p(550). Their masses do not vanish in the chiral 
limit 0)- our range of temperature, they behave like a dilute gas. They 

are of course exponentially suppressed, but their effects become more and more 
important when the temperature is increased. In the order parameter they generate 
a contribution of 0.5% at T = 100 MeV with respect to that of the pions, whereas 
at T = 150 MeV it becomes of the order of 10%. At T ~ 160 MeV, the mean 
distance between the massive states is approximately 1.6 fm and the number of 
massive particles per unit volume is the same as the numbers of pions. This gives a 
limitation to the GhPT approach. In the present calculation this limitation shows up 
at a similar temperature: the square of the velocity of propagation of the Goldstone 
bosons turns negative. Our results are meaningless beyond that point. 
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A Properties of the T-dependent integrals 


A direct computation of the closed-eye,the guimbard and the sunset graphs brings 
different complicated integrals. But taking into account the symmetry of the inte- 
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grands one can reduces them to two independent integral forms. A short notation 
is hrst introduced: 

{fiQ, h, h)) ■= i J ^ ^ 27 !) 2 d^ - h - fca) f{q, h, ^ 2 ). (A.l) 

Because (/(g, fci, ^ 2 )) = (/(g, ^ 2 ,^ 1 )) = (/(g, g - - ^ 2 , ^ 2 )) and that (kj) = 

M^I{q,T) — A^(0), the expressions for the closed-eye, the guimbard and the sunset 
can be simplihed into ( p.22| , |3.23i|3.24|) . 

These two integrals /(g, T) and I/i^iq, T) diverge in d = 4 dimensions. They have 
to be regularized. They contain three types of contributions. The T-independent 
terms are not of interest here. The parts involving one Bose distribution diverge for 
d = 4 and the ones with a product of two statistical weights are hnite for d = 4. 


The hrst function of interest can be written as: 


/(g, T) = J(g, T) - 6 Wm(T) - ^ 


where A/m(T), already introduced in ( |3.13|) , is 

d^k 


(2 

and I{q,T) is a hnite integral given by 

d^k 


r d^k 

Mm{T) = j T^5{k^ - M^) ubM, 


I{q,T) = 3j 
+3 
+i 


(27r)3 

r d^k 


5{k'^ - M‘^)nB{oOk)J{{q + kf) 
dik'^ - M^)nB{uJk)K{q + k,T) 


(27r 

- M‘^)nB{u}ki)d{kl - M‘^)nB{uJk2) 

^((g -h- k2f - M^)nB{ujq-ki-k2] 


(A.2) 


(A.3) 


(A.4) 


The second function we are interested in is 
V(g,^) = hu{q,T) 

+A {Nm(T) \g,u(\q^ - \m'^) - - Y N^,{M, T)) (A.5) 

A^M(r) [4 + M^) - J ^ T)) , 


Ibvr^ 

where is 


c k 

iV^,(M,T) = J -^k^kJik^ - M^)nBiuJk), 


(A.6) 
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and the finite integral T) is 

r 

hM.T) = j ^^5{e-M^)nBMJ{{q + kf) 

f 5 2 2 

+ ^q^,qu + g(g/.A;^ + q^k^) 

2M^ (gAt + fc^)(g^ + fc^) ] 

^ 2 6 3 ; 3 (g + A;)2 J 

r d^k 

+ j {-2k^Mq + fc, T) 

+/t'(g + k,T)[Ak^ku + q^qu + 2(g^fci. + g^/c^)]} (A.7) 

+* J ~ M‘^)nB{uJk^)5{kl - M‘^)nB{uJk2) 

^((<? - ^1 - hf - M‘^)nB{ujg-ki-k2) 

{Aki^ki^ + g^g^ - 2{q^ki^ + q^ki^) + 2ki^k2y} . 


In (|A.4|JA.7| ), three one-loop functions have been introduced to get a simpler 
representation. The usual T = 0 one 


W) = 


1 


Ihvr^ 


dx\n{l — g^x(l — x)/M‘^), 


together with temperature dependent ones: 

d^k 


K{q,T) = z 


( 2 vr: 


^6{k‘^ - M^) nB{uJk) A(g - k) 


1 


rOO 

/ dk —nB(a;fc)ln 
Jo OJh 


IhTT^I^ 


k 


{q^-2kW-Aixiqi 
{q^ + 2k\^Y - Aul ql 


and 


KM^T) = i 


SttI^ 
f d^k 


/ dk — UBiuJk) 
Jo OOk 


(27r) 


■6{k‘^ - M^) UBiuJk) A(g - k)k^ 


(A.8) 


(A.9) 


= ^(n,[^AfMiT) + ^K{q,T)-q^Koiq,T)] 

+ ^K{q,T) - qoKo{q,T)]^ , 

where n = (1, 0, 0, 0). An integral representation for KQ{q, T) can be given: 


(A. 10) 


Ko{q,T)^ 


Ibvr^l^ Jo 


dk k UBioJk) In 


g^ -4(o;fc go - k\^f 
g4 - 4(a;fc go + k\^Y 


(A.ll) 
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Note that because of the symmetry of the problem under spatial rotations: 


and 


N,k{M,T) = 1 5,k{Noo{M,T) - M^Mm{T)), 

(A.12) 

IoM,T) = qiR{q,T)lA 

(A. 13) 

Iik{q,T) = Sik Ii{q,T) + qiqk l2{q,T). 

(A.14) 


The examined observables are dehned on-shell at ^ = 0 and = M. In that 
special case, a simplihed representation can be given. Note that because Uq > M, 
the different sums of T-independent one-loop functions appearing in the integrands 
can be rewritten as: 


J{2M^ + 2Muk) + J{2M^ - 2Muk) 
J{2M^ + 2MuJk) - J{2M^ - 2Muk) 


1 

1 


(2 + ^ \n(j{k) 

\ \k\ 

[2 + — \naik) 

\ \k\ 


where 


a{k) := 


\/ uik + M — ■\/ ujk — M 
\/Uk + M (Jk — M 


The following representations can be obtained: 

2 /• ^ r. 2 

\qo=M,q=0 


Re[J(g,T)] 


(2+i^log<rW 


2 poo pi 

+ T 7 r^ / dk da - nsioo) nsiuJa) In 

lOTT^ Jo Jo 




(A.15) 


(A. 16) 


(A. 17) 

1 “h O; 

1 — a 


R-e[Joo(q', T)] 


qQ=M,q=0 


1 poo f 

/ dk — nB{uJk)\{5k‘^+ 7M‘^) (A.18) 

9 o 7 r^ Jo Uk [ 

(2 + ‘7^ loga{k)^ + 2uJkM ( 2 - ^ log a(A;)^ | 

poo pi 

+ / dk da ak^nBiuJk) nsiuJa) Inrfo;,/c). 

Idtt^ Jo Jo 


The expressions given above involve some new functions: 


oJa '■= V a^k"^ + 


(A.19) 


23 



r(a, k) := 


2k^a^ + M^(l + a^) — 2auJkOOc 


(A.20) 


2kP‘a^ + M2(l + a^) + 2au}k^a 
In general these integrals cannot be algebraically evaluated, even if the mass is zero. 


Finally Iik{q,T) 
because of (|A.14|) : 


qo=M,q=0 


can be expressed in term of other known functions 


hk{q,T) 


qo=M,q=0 


\ Sik {loo{q,T)-M^I{q,T) +AA^(T)) 


qo=M,q=0 


(A.21) 


B Some functions in the chiral limit 


As already mentionned all our integrals are in fact functions of M/T. It is however 
not so easy to evaluate the desired ones in the chiral limit. Some clever tricks can 
be found in @]. The value of the four functions (|3.13|JA.3|JA.17||AT^) in the chiral 
limit M T is: 


Mm{T) 


T2 


m=0 


2'k‘^ Jo 
MT 

12 ~ ^ 


dt t nB{tT) 


M 


27r2 T 


/ dt — 1= 
n t \/W 


= + 0 (^ln—) 




(B.l) 


Noo{M,T) 


m=0 


2^4 


poo 

/ dt t^ usitT) + 0{ 
Jo 


T' 


Re[J(g,T)] 


Re[Joo(g,T)] 


5=0,771=0 


30 + ( 7 - 2 )’ 


(B.2) 

M f 

f)), 

(B.3) 





771=0 

(B.4) 


.oo 


1 + q ;2 


H- - I dt da ar nsitT) UBiatT) In ■ 

Stt^ Jo Jo a^k^ 
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1 


= T 


/ 1 M 

ll44‘"Y 


-In 2 H- 

144 1728 


0.0073). 0(£). 


The numbers appearing in the previous expressions are due to numerically eval¬ 
uated integrals that come directly from the representations shown in Appendix A. 


K{q,T) 


g=0,ni=0 


6M 




144Mir2 






(B.5) 


_d_ 

dqo 


Re[^oo(9t r)] 


q=0,7h=0 


_ 


Svr'^ 


2M 


C(3) 


(A/m(T)L^o) 

3T3 


48M7r2 


Noo{M,T)\^^, 


/ dtf nsitT) 

lOTT^ JO 

1 \ M, 

+ 0(—In—), 


487r 


T T 


(B.6) 


The cancellation of the 1/M terms in the last integral in the chiral limit is very 
important for the safe of the GOR relation at hnite temperature and the term 
is responsible of the corrections (|6.16|). 
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Figure Captions 


Fig. 1 : The Keldysh path of integration in the complex t plane. The arrows indicate 
the time ordering on each part of the contour. Only Ci and C 2 are relevant in the 
generating functional. 


Fig. 2 : The Feynman diagrams necessary to compute the axial two-point Green’s 
function to two loops. The various vertices correspond to the part of the effective 
Lagragian ( p.2|) involved: the ’’dot” for the vertices, the ”4” (resp. ”6”) for 
those of (resp. The wiggled lines represent an external held, whereas the 

plain lines are the thermal propagators. The various parts of the RTF propagator 
and the crossed graphs are not explicitly given. 


Fig. 3 : The effective pion mass at non-zero temperature. The dashed-doted curve 
represents the trivial result at the tree level, the dashed one the one-loop computa¬ 
tion and the full one the two-loop approximation. 
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Fig. 4: The real part of the ’’temporal” pion decay constant at non-zero temper¬ 
ature. The dashed-doted curve represents the trivial result at the tree level, the 
dashed one the one-loop computation and the full one the two-loop approximation. 


Fig. 5 : The T-dependent pion mass to two loops in the chiral limit (dashed curve) 
and in the physical case (full one). 


Fig. 6; The real part of the ’’temporal” pion decay constant to two loops in the 
chiral limit (dashed curve) and in the physical case (full one). 


Fig. 7: The difference between the real parts of the ’’temporal” and ’’spatial” pion 
decay constants to two loops in the chiral limit (dashed curve) and in the physical 
case (full one). In the chiral limit this quantity is related to the speed of the pions. 
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